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Abstract 

In this paper we present a decision procedure for computing 
order linear ODEs, that is, solutions for the classes of hypergeometric equations constructed from the 
3F2, 2F2, 1F2 and 0F2 standard equations using transformations of the form x —* F(x), y — > P(x)y, 
where F(x) is rational in x and P(x) is arbitrary. A computer algebra implementation of this work is 
present in Maple 12. 

Introduction 

Given a third order linear ODE 

y'" + c 2 y" + c iy ' + c y = (1) 
where y = y{x) is the dependent variable and the Cj = Cj(x) are any functions of x such that the quantities 1 

~ 3~ Io = — -^T + — - Co (2) 

are rational functions of x, the problem under consideration is that of systematically computing solutions for 
(1) even when no Liouvillian solutions exist 2 . Recalling, Liouvillian solutions can be computed systematically 
[1] and implementations of the related algorithm exist in various computer algebra systems. The linear 
ODEs involved in mathematical physics formulations, however, frequently admit only non-Liouvillian special 
function solutions, and for this case the existing algorithms cover a rather restricted portion of the problem. 

The special functions associated with linear ODEs frequently happen to be particular cases of some gen- 
eralized hypergeometric p F q functions [2] . One natural approach is thus to directly search for p F q solutions 
instead of special function solutions of one or another kind, and this is the approach discussed here. Related 
computer algebra routines were implemented in 2007 and are now at the root of the Maple (release 12) [3] 
ability for solving non-trivial 3rd order linear ODE problems. 

The approach used consists of resolving an equivalence problem between a given equation of the form 
(1) and the four standard p F q differential equations associated to third order linear ODEs, that is, the 3F2, 
2F2, 1F2 and 0F2 equations [4], respectively: 



x 7i and Io are invariant under transformations of the dependent variable of the form y(x) — > P(x) y(x), P arbitrary. 
2 Expressions that can be expressed in terms of exponentials, integrals and algebraic functions, are called Liouvillian. The 
typical example is exp( J R(x), dx) where -R(x) is rational or an algebraic function representing the roots of a polynomial. 
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where {a, /?, 7, <5, 77} represent arbitrary expressions constant with respect to x. The equivalence classes are 
constructed by applying to these equations the general transformation 3 

x^F(x), y^P(x)y (4) 

where P(x) is arbitrary, with the only restriction that F(x) is rational in x, resulting in rather general ODE 
families. When the equation being solved belongs to this class, apart from providing the values of F(x) and 
P{x) that resolve the problem, the algorithm systematically returns the values of the (five, four, three or 
two) p F q parameters entering each of the three independent solutions. 

It is important to note that the idea of seeking hypergeometric function solutions for linear ODEs or 
using an equivalence approach for that purpose is not new, although in most cases the approaches presented 
only handle second order linear equations [6, 7, 8, 9]. An exception to that situation are the algorithms 
[10, 11] for computing p F q solutions for third and higher order linear ODEs, and a similar one implemented 
in Mathematica [12]. It is our understanding, however, that the transformations defining the classes of 
equivalence that those algorithms can handle are restricted to x — > ax b , y — > P(x)y, with a and b constants, 
not having the generality of (4) with rational F(x) presented here. 

Apart from concretely expanding the ability to solve third order linear ODEs, the decision procedure 
being presented generalizes previous work in that: 

1 . The ideas presented in [9] , useful for decomposing two sets of invariants into each other, were extended 
for third order equations and elaborated further. 

2. The classification ideas presented in [9] for second order linear equations were extended for third order. 

3. When the p F q parameters are such that less than three independent p F q solutions exist, instead of 
introducing integrals [11], MeijerG functions are used to express the missing independent solutions. 

The combination of items 1 and 2 resulted in the new ability to solve the p F q ODE classes generated by 
transformations as general as (4) with F(x) rational. Item 3 is not new 4 , though we are not aware of literature 
presenting the related problem and solution. Altogether, these ideas and its related algorithm permit the 
systematic computation of three independent solutions for a large set of third order linear equations that we 
didn't know how to solve before. 



1 Computing hypergeometric solutions 

To compute p F q solutions to (1) the idea is to formulate an equivalence approach to the underlying hyper- 
geometric differential equations, that is, to determine whether a given linear ODE can be obtained from one 

3 The problem of equivalence under transformations {x — ► F(x), y — ► P(x) y + Q(x)} for linear ODEs can always be mapped 
into one with Q(x) = 0, sec [5]. 

4 Mathematica 6 also uses MeijerG functions as described in item 3. 
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of the p Fg ODEs (3) by means of a transformation of a certain type. If so, the solution to the given ODE is 
obtained by applying the same transformation to the solution of the corresponding p F q equation. 

The approach also requires determining the values of the hypergeometric parameters {a, /?, 7, <5, 77} for 
which the equivalence exists, and it is clear that the bottleneck in this approach is the generality of the 
class of transformations to be considered. For instance, one can verify that for linear transformations of the 
form (4) with arbitrary F(x), in the case of second order linear ODEs, the problem is too general in that 
the determination of F(x) requires solving the given ODE itself [13], making the approach of no practical 
use. This has to do with the fact that in the second order case, any linear ODE can be obtained from any 
other one through a transformation of the form (4). The situation for third order equations is different: the 
transformation (4) is not enough to map any equation into any other one 5 [14], so that its determination 
when the equivalence exists is in principle possible. By restricting the form of F{x) entering (4) to be rational 
in x the problem becomes tractable by using a two step strategy: 

1. Compute a rational transformation R(x) mapping the normal form of the given equation 6 into one 
having invariants with minimal degrees (defined in sec. 3). 

2. Resolve an equivalence problem between this equation with minimal degrees and the standard p F q 
equations (3) under transformations of the form discussed in [9], that is 

{ax k +b) 

X — > 7 ; --, V^PlXjV 5 

with P(x) arbitrary and {a, b, c, d, k} constants with respect to x. In doing so, determine also the 
parameters {a, (3, 7, <5, r]} of the p F q or MeijerG functions entering the three independent solutions. 

The key observation in this "two steps" approach is that a transformation of the form (4) with rational F(x) 
mapping into the p F q equations (3), when it exists, it can always be expressed as the composition of two 
transformations, each one related to each of the two steps above (see sec. 3), because (3) have invariants with 
minimal degrees. The advantage of splitting the problem in this way is that the determination of R(x) in 
step one, and of the (up to five) p F q parameters in step two, as well as of the values of {a, b, c, d, k} entering 
(5), is systematic (see sec. 2 and sec. 3), even when the problem is nonlinear in many variables. 

2 Equivalence under x — > (a x k + b)/(cx k + d), y — > P(x) y 

This type of equivalence is discussed in [9] and generalized here for third order ODEs. Recalling the main 
points, these transformations, which do not form a group in the strict sense, can be obtained by sequentially 
composing three different transformations, each of which does constitute a group. The sequence starts with 
linear fractional - also called Mobius - transformations 

— 1 j, 

(6) 



cx + d' 

is followed by power transformations 

x->x k , (7) 
and ends with linear homogeneous transformations of the dependent variable 

y^Py. (8) 



5 Thereforc there exist enough absolute invariants to formulate the equivalence problem under (4) - see sec. 3. 

6 The coefficients of y' and y in the normalized equation are the invariants Ii and Iq defined in (2), assumed to be rational. 
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2.1 Equivalence under transformations of the dependent variable y — > P(x)y 

Transformations of the form (8) can easily be factored out of the problem: if two equations of the form (1) 
can be obtained from each other by means of (8), the transformation relating them is computable directly 
from these coefficients. For that purpose first rewrite both equations in normal form using 

— / C2(x)/3 dx / n \ 

V -» ye J (9) 
and the transformation relating the two hypothetical ODEs - say with coefficients Cj and Cfc, when it exists, 
is given by y - y e /(«*(*)-&(*))/3«fa 

2.2 Equivalence under Mobius transformations, singularities and classification 

Mobius transformations preserve the structure of the singularities of (1). For example, all of the oF 2 , iF 2 and 
2F2 hypergeometric equations in (3) have one regular singularity at the origin and one irregular singularity 
at infinity, and after transforming them using the Mobius transformations (6), they continue having one 
regular singularity and one irregular singularity, now respectively located at 7 — b/a and —d/c. 

In the case of the 3F2 differential equation (the first listed in (3)), under (6) the three regular singularities 
move from {0, 1, 00} to {— b/a, —d/c, (d—b)/(a — c)}. So from the singularities of an ODE, not only one can 
tell with respect to which of the four differential equations (3) could the equivalence under (6) be resolved, 
but also one can extract the values of the parameters {a, b, c, d} entering the transformation (6). 

More generally, through Mobius transformations one can formulate a classification of singularities of the 
linear ODEs "equivalent" to the third order p F q equations (3) as done in [9] for second order p F q equations. 
So, for each p F q family obtained from (3) using (6), a classification table can be constructed based only on: 

• the degrees of the numerators and denominators of the invariants (2); 

• the presence of roots with multiplicity in the denominators; 

• the possible cancellation of factors between the numerator and denominator of each invariant. 

With this classification in hands, from the knowledge of the degrees with respect to x of the numerator 
and denominator of the invariants (2) of a given third order linear ODE, one can determine systematically 
whether or not the equation could be obtained from the 3F2, 2F2, iF 2 or 0F2 equations (3) using (6). 

2.3 Transformations x — > F(x) and equivalence under x — > x k 

Changing x — ► F(x) in (1), the new invariants Ij can be expressed in terms of the invariants (2) of (1) by 

h(x) = F' 2 h(F)-2S(F) 

(10) 

Io(x) = F'F"h(F) + F' 3 I Q (F) - S(F)' 

where S(F) is the Schwarzian [15] 

F'" 3 f F" s 2 



The form of S(F) is particularly simple when F{x) is a Mobius transformation, in which case S(F) = 0. 

Regarding power transformations F(x) = x k , unlike Mobius transformations, they do not preserve the 
structure of singularities; the Schwarzian (11) is: 

s(* k ) - y^. (12) 

From (10) and (12), for instance the transformation rule for h(x) becomes 



7 When either a or c are equal to zero, the corresponding singularity is located at 00 
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x 2 h{x) + l = k 2 {{x k fh{x k ) + 1) . (13) 
Generalizing to third order the presentation of shifted invariants in [9] , we define here 

Ji(x) = x 2 h(x) + 1, 

(14) 

J 2 (x) = x 3 I (x) + x 2 I 1 (x). 
From (10) rewritten in terms of these J n (x), their transformation rule under x — ► x k is given by 

Ji(x) = fc 2 Ji(x fe ), 

(15) 

J 2 {x) = k 3 J 2 (x k ). 

The equivalence of two linear ODEs A and B under x — > x k can then be formulated as follows: Given 
the shifted invariants J n ,A(x) and J,i,b(x), computed using their definition (14) in terms of I n {x) defined 
in (2), compute fc^ and ks entering (15) such that the degrees of J u ,a{x) and J n ,B{x) arc minimal. From 
the knowledge of x —* x kA and x — > x kB , respectively leading to J Uy A and J Uy B with minimized degrees, 
equations A and B are related through power transformations only when J n ^A = J n ,B and, if so, the mapping 
relating A and B is x — ► x kA ~ kn . Finally, the computation of k simultaneously minimizing the degrees of 
the two J n {x) in (15) is performed as explained in section 3 of [9]. 

3 Mapping into equations having invariants with minimal degrees 

The decision procedure presented in the previous section serves for systematically solving well defined families 
of pF q 3rd order equations for which no solving algorithm was available before to the best of our knowledge. 
However, the restriction in the form of F{x) entering (4) to the composition of Mobius with power transfor- 
mations is unsatisfactory: for linear equations of order higher than two, (4) does not map any linear equation 
into any other one of the same order and so the problem is already restricted 8 . 

As shown in what follows, one possible extension of the algorithm is thus to consider the general trans- 
formations (4) restricting F{x) to be a rational function of x. For that purpose, instead of working with 
invariants Ij under y — > P(x) y we introduce absolute invariants Lj under {x — > F(x), y — > P(x) y}: 

(6rr" + 9hr 2 - 7r' 2 ) 3 T (27/i'r 3 - I8hr 2 r' + 56r' 3 - 72r"r'r + 18r"V 2 ) 
L i = ~8 > ^ = - 4 ; (16) 

where r = I\ — 2I is a relative invariant of weight 3 [16]. Under (4), Li transforms as Li(x) — ► Li(F{x)) and 
(16) can be inverted using as intermediate variables the relative invariants s — (L 2 L\) / L\ , and t — Li/s 3 : 

st 3 - 6t"t + It' 2 T (s' - 9)t 4 + t'st 3 - 6t"'t 2 + 20t"t't - Ut' 3 

h = m 2 ' /o ~ 18^ ' (17) 

Thus, any canonical form for the Li that can be achieved using (4) automatically implies on a canonical 
form for the h and so for the ODE (1). The canonical form we propose here is one where the Li have 
minimal degrees, that is, where the maximum of the degrees of the numerator and denominator in each of 
the Li(x) is the minimal one that can be obtained using a rational transformation x — ► F{x). This canonical 
form is not unique in that it is still possible to perform a Mobius transformation (6), that changes the Li 
but not their degrees. 

The equivalence of two linear ODEs A and B under (4) with rational F{x) can then be formulated 
by rewriting both equations in this canonical form, where the invariants L, of each equation have minimal 
degrees, followed by determining whether these canonical forms are related through a Mobius transformation. 



8 The equivalence problem for linear equations of order n involves a system of n — 1 equations and invariants Ij(x), that 
includes the equivalence function F(x). When n > 2, eliminating F(x) from the problem results in an interrelation between 
the Ij so that the equivalence is only possible when these relationships between the Ij hold [14] . 
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In the framework of this paper, B is one of the hypergeometric equations (3), all of them already in 
canonical form in that the corresponding Li already have minimal degrees. Hence, the equivalence of A, of 
the form (1), and any of B of the form (3) requires determining only a canonical form for A (the rational 
function F(x) minimizing the degrees of the Li of A), followed by resolving an equivalence under Mobius 
transformations between this canonical form and any of the equations (3), done as explained in sec. 2.2. 

The key computation in this formulation of the equivalence problem under (4) is thus the computation 
of a rational F(x) that minimizes the degrees of the Li of A. The computation of F(x) can clearly be 
formulated as a rational function decomposition problem subject to constraints: "given two rational functions 
Li(x),i = 1..2, find rational functions Li{x) and F(x) satisfying Li — Li o F and such that the rational 
degree of F is maximized" (and therefore the degrees of the canonical invariants Li are minimized). In turn, 
this type of function decomposition associated to "minimizing the degrees" of the Li can be interpreted 
as the reparametrization, in terms of polynomials of lower degree, of a rational curve that is improperly 
parameterized, as discussed in [17], where an algorithm to perform this reparametrization is presented. 

One key feature of the algorithm presented in [17] is that it reduces the computation of F(x) to a sequence 
of univariate GCD computations, avoiding the expensive computation of bivariate GCD. However, it is not 
clear for us whether the prescriptions in [17] (at page 71) for mapping the bivariate GCDs into univariate 
ones is complete. We also failed in obtaining a copy of the computer algebra packages FRAC [18] or Cadecom 
[19] that contain an implementation of the algorithm presented in [17]. Mainly for these reasons, and without 
the intention of being original, we describe here a slightly modified version of the algorithm presented in [17]. 

3.1 An algorithm for computing x — > F(x) minimizing the degrees of the Lj(x) 

Let Li(x) — Li(F(x)) — Ni(x)/ Di(x),i — l..n, and F(x) — p(x)/q(x), where the Li have minimal degrees, 
Ni is relatively prime to Dj and p is relatively prime to q. Construct polynomials 

Qi(x,t) = numerator (L< (x) - L t (t)) = Ni(x)Di(t) - Ni(t)Di(x), (18) 
and let P(x,t) be the bivariate GCD of these Qi(x,t). Consequently 

P(x, t) = numerator(F(a;) - F{t)) = ^ P l {x) t l = p(x) q(t) - p(t) q(x), (19) 

i 

The coefficient Pi(x) of each power of t in P(x,t) is a linear combination of p(x) and q(x), and because 
the quotient of any two relatively prime of these linear combinations is fractional linear in F(x), so is the 
quotient of any two relatively prime Pi{x). Finally, because F(x) is defined up to a Mobius transformation 
we can take that quotient itself - say, Pi(x)/Pj(x) - as the solution F(x). 

The slowest step of this algorithm is the computation of the bivariate GCD between the Qi(x,t) that 
determines the function P(x, t) from which the Pi(x) are computed. It is possible however to avoid computing 
that bivariate GCD, using a small number of univariate GCD computations instead. 

For that purpose, notice first that what is relevant in the Pi(x) is that they are linear combinations of p(x) 
and q(x). Now, we can also obtain linear combinations of p{x) and q{x) by directly substituting numerical 
values tk for t into P(x,t), and from there compute F(x) as the quotient, e.g., of P(x,to)/ P{x,t\). The key 
observation here is that these P{x,tk) can also be obtained by substituting t = tk directly into the Qi(x,t) 
followed by computing the univariate GCD of Qi(x, tk) and Qi{x, tk) 9 , avoiding in this way the computation 
of the expensive bivariate GCD leading to P(x,t). 

Repeating this process with another i-value gives a second, in general different, such linear combination 
of p(x) and q(x), with F being the resulting quotient of two of these linear combinations obtained using 
different values of t. The rest of the algorithm entails avoiding invalid i-values at the time of substituting 
t = tk and this is accomplished by considering different tk until the following conditions are both satisfied: 

9 For example, suppose the ^-solutions of P(x,t) = are x = Xj(t),j = l..m, i.e., P(x,t) = -P m (t)Y[ .x-Xj(t). Then 
each x = Xj(to) is a solution of both Qi(x,to) = and Q2(^>to) = 0- F° r most values of to (all but a finite set in fact) these 
Xj(to) will be the only such common solutions, and therefore the GCD of Qi(x,to) and Q2(x,to) is in fact P(x,to)- 
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1. The two P(x,to), P(x,t\) whose quotient gives the solution F(x) must be relatively prime. 

2. The degree of F must divide the degrees of each Li, i = l..n. 

4 Summary of the p F q approach for third order linear ODEs 

The idea consists of assuming that the given linear ODE is one of p F q equations (3) transformed using (4) 
for some F(x) rational in x and P(x) arbitrary and for some values of the pFq parameters. Resolving the 
equivalence is about determining the F(x), P{x) and the values of the p F q parameters {a, (3, 7, 5, 77} such 
that the equivalence exists. An itemized description of the decision procedure to resolve this equivalence, 
following the presentation the previous sections, is as follows. 

1. Rewrite the given equation (1) we want to solve, in normal form 

y"' = h(x)y' + I (x)y (20) 
where the invariants I n (x) are constructed using the formulas (2). 

2. Verify whether an equivalence of the form {x — > (ax k + b)/(cx k + d), y — > P{x) y} exists: 

(a) Compute J n (x), the shifted invariants (14), and use transformations x — ► x k to reduce to the 
integer minimal values the powers entering the numerator and denominator; i.e., compute k and 
J n (x) in (15). 

(b) Determine the singularities of the J n (x) and use the classification of singularities mentioned in 
section 2 to tell whether an equivalence under Mobius transformations to any of the 3 F 2 , 2F2, 
1F2 or F 2 equations (3) exists. 

(c) When the equivalence exists, from the singularities of the two J n (x) compute the parameters 
{a, b, c, d} entering the Mobius transformation (6) as well as the hypergeometric parameters 
{a, (3, 7, 5, 77} entering the p F q equation (3). 

(d) Compose the three transformations to obtain one of the form 

ax k +(3 
x^ = -, y—>P(x)y 

mapping the p F q equation involved into the ODE being solved. 

3. When the equivalence of the previous step does not exist, perform step 1 in the itemization of section 1, 
that is, compute the absolute invariants (16) and compute a rational transformation R(x) minimizing 
the degrees of the invariants (16) of the given equation 

(a) When R(x) is not of Mobius form, change x — > R(x) rewriting the given equation in canonical form 
and re-enter step (2) with it, to resolve the remaining Mobius transformation and determining 
the values of the p F q parameters. 

4. When either of the equivalences considered in steps (2) or (3) exist, compose all the transformations 
used and apply the composition to the known solution of the p F q equation to which the equivalence 
was resolved, obtaining the solution to the given ODE. 
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5 Special cases and MeijerG functions 



Giving a look at the series expansion of any of the 3F2, 2F2, 1F2 or 0F2 functions one can see that there 
are some different situations that require special attention at the time of constructing the three independent 
solutions to (1). Consider for instance the standard 0F2 equation and its three independent solutions, 



x x z x z (21) 

y = F 2 ( ; a,0; x) d + x l ~^ F 2 { ; 2 - /?, 1 + a - (3; x) C 2 + x 1 - Q F 2 ( ; 2 - a, 1 - a + (3; x) C 3 

where the d are arbitrary constants. Expanding in series the first 0F2 function entering this solution we get 

1 + + 2af3(a + l)(l+pf 2 + &a(3 (a + 1) (1 + (3) (a + 2) {(3 + 2)^ + ° ^ (22) 

This series does not exist when a or (3 are zero or negative integers, and the same happens when the 
pF q parameters entering any of the other two independent solutions is a non-positive integer. By inspection, 
however, one of the three p F q functions entering the solution in (21) always exists, because there are no a 
and (3 such that the three 0F2 functions simultaneously contain non-positive integer parameters. 

Consider now the second independent solution, x 1 ~ /3 F 2 ( ; 2 — (3, 1 + a — /?; x): when (3 = 1 it becomes 
equal to the first one and so we have only two independent p F q solutions. In the same way, when a = 1 
the first and third solutions entering (21) are the same and when a = (3 the second and third solutions are 
the same. And when the two conditions hold, that is a = (3 = 1, actually the three solutions are the same. 
Notwithstanding, in these cases too one of the three 0F2 solutions always exists. 

The same two type of special cases exist for the 1F2, 2F2 and 3F2 function solutions and the problem 
at hand consists of having a way to represent the three independent solutions to (1) without introducing 
integrals or iterating reductions of order 10 . For this purpose, we use a set of 3 MeijerG functions for each 
of the four p F q families that can be used to replace the missing p F q solutions in these special cases. The 
key observation is that at these special values of the last two parameters of the p F q functions the MeijerG 
replacements exist, satisfy the same differential equation and are independent of the available p F q function 
solutions. A table with these 3 x 4 = 12 MeijerG function replacements is as follows: 



Table 1: MeijerG alternative solutions to the p F q equations 



p F q family 



MeijerG functions 



oF 2 ( ; a, 13; x) 
iF 2 (a; (3,r, x) 
2 F 2 (a,/3; (5,7; x) 
3 F 2 (a,/3,7; S,rj; x) 



^3 ( •' ■ 



0,l-a,l-/3 
1-a 

0,1-/3,1-7 

1-/3, 1-a 
0,1-7,1-5 



r 2,3 f l-/3,l-a,l- 7 \ 
^3,3 \ x > I 0,1-5,1-^ J 



G\ 



.3.0 
0.3 



G 3 {1 ( 



G 



3.2 
2,3 



0,l-a,l-/3 
1-q 

0,1-7,1-/3 

1-/3, 1-a 
0,1-7,1-5 



^3,3 ( l-/3,l-a,l- 7 \ 
U 3,3 ^ x i I 0,1-5,1-)) ) 



U 1.3 



[ x ' 



1-a, 1-/3,0 
1-a 

1-7,1-/3,0 

1-/3, 1-a 
1-7,1-5,0 



r 2,3 ( l-/3,l-a,l- 7 \ 
U 3,3 \ x > I 1-5,1-?;, J 



6 Examples 



Equivalence under power composed with Mobius transformations for the F 2 class 

Consider the third order linear ODE 

10 Recall that given two independent solutions, it is always possible to write the third one in terms of integrals constructed 
with the two existing solutions, and in the case of a single solution it is still possible to reduce the order to a second order linear 
equation that may or not be solvable. 
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(37 + 2^ + 6,-108^) 

V ~ \2x(x + l)(x-l) V (26) 

(2 (, + 6) (11/2- n) + (36^ + 294+ 12 M ).t 2 - 360x 4 ) , 16 
+ 24 s 2 (x + l) 2 {x-lf V x (x + l) 4 (x - l) 4 V 

This equation has two regular singularities at {0, oo} and two irregular singularities at { — 1, 1}. Following 
the steps mentioned in the Summary, we rewrite the equation in normal form and, in step 2. (a), compute 
the value of k leading to an equation with minimal degrees entering J n {x) in (15). The value of k found is 
k = 2 so the equation from which (23) is derived changing x — > x 2 is 

_ (6, + 2 M + 73- 144a;) 

V ~ 2Ax(x-l) V (24j 

(2 (i/ + 8)(/x+l/2) - (48z/+16^ + 584)x + 576a; 2 ) , 2 

96 a; 2 {x-lf V ~ x 2 (x - l) 4 V 

and has invariants with minimal degrees with respect to power transformations. In step 2.(b), analyzing 
the structure of singularities of (24) we find one regular singularity at the origin and one irregular at oo. 
Using the classification discussed in section 3.2 based on the degrees with respect to x of the numerators 
and denominators of the invariants of (24) as well as the factors entering these denominators the equation 
is identified as equivalent to the 0F2 class under Mobius transformations (6). So we proceed with step 2.(c), 
constructing the Mobius transformation and computing the values of the hypergeometric parameters {/i, 1/} 
entering the 0F2 equation in (3) such that the equivalence under Mobius exists, obtaining: 

a = i//4 + 2, /3 = /i/12+ 1/24, M := x -> (25) 

x — 1 

Composing M above with the power transformation used to obtain (24) and using the values above for a 
and P, in step 4 we obtain the solution of (23) 

y(x) = F 2 ( ;i//4 + 2,^/12 + 1/24; 2 ^—j)Ci 

+x-( 2 +^ 2 ) ( x 2 _ lf +v/i) Q F 2 ( ; -iz/4,^/12- 23/24; ^-y)C 2 (26) 

+a .(23/12-„/6) ( x 2 _ ^ (M/12-23/24) ^ . ^ ^ ^ ^ + ^ _2x^_ } ^ 

Meijerg functions and equivalence under rational transformations for the iF 2 class 

Consider the following equation, with no symbolic parameters and only integer powers 

/// (6 +12 a; -15 a; 2 -6 a; 3 ) „ 

V ~ x(l + x-x 2 )(x + 2) V (2?) 

(16 + 48 a; + 36 a; 2 - 20 a; 3 + 9 a; 4 + 81 a; 5 - 20 a; 6 - 30 a; 7 -6 a; 8 ) , 
+ a; 4 (a; + 2) 2 (l + x-a; 2 ) 2 V 



(1 + x — x 2 ) x 



5 
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Following steps 1 and 2 in the Summary, we confirm that there exists no equivalence under (5), so in step 3 
we search for a rational transformation minimizing the degrees of the invariants (16), finding 



R(x) = x 2 /(l + x) 
Therefore (27) can be obtained by changing variables x — > R{x) in 



y = — 



9x + 6x 2 ) „ (l 

■ y" + — 



2x + 6x 2 -6x 3 ) 



1 



y - 



■y 



(28) 



(29) 



x(x — 1) x 3 (x — 1) {x — iy x" 1 

This equation 11 thus has invariants with minimal degrees, and has one regular singularity at 1 and one 
irregular at the origin. According to the classification in terms of singularities (29) admits an equivalence 
under Mobius transformations to the p F q equations (1F2 case) and hence is solved in the iteration step 3. (a) 
mentioned in the summary. When constructing the p F q solutions to (29), however, we find that the iF 2 
parameters in the second list are both equal to 1, so only one 1F2 solution is available, and hence two of the 
MeijerG alternative solutions presented in the table (5) are necessary, resulting in 



y - o? i{ , 1, — — 2 — ) °i + G o,2 



0,0 ) C2 + gI'I 



1 



0,0,0 



(30) 



Note that the first p F q function is a Fi. This is due to the automatic simplification of order that happens 
when identical parameters are present in both lists of a iF 2 function; this 0F1 can also be expressed in terms 
of Bessel functions. 



Conclusions 

In this work we presented a decision procedure for third order linear ODEs for computing three independent 
solutions even when they are not Liouvillian or when the hypergeometric parameters involved are such that 
only two or one p F q solution around the origin exists. This algorithm solves complete ODE families we didn't 
know how to solve before. 

The strategy used is that of resolving an equivalence problem to the 3F2, 2F2, iF 2 and 0F2 equations, 
and in doing so, two important generalizations of the algorithm presented in [9] were developed. First, 
the classification according to singularities and the use of power composed with Mobius transformations, 
presented in [9] for 2nd order equations, was generalized for third order ones. Second, the idea of resolving the 
equivalence mapping into an equation with invariants with "minimal degrees under power transformations" 
was generalized by determining a transformation mapping into an equation having invariants with "minimal 
degrees under general rational transformations" . This permits resolving a much larger class of p F q equations, 
defined by changing variables in (3) using {x — > R(x), y — ► P(x) y} where R(x) is a rational function. 
Symbolic computation routines implementing this algorithm were integrated into the Maple system in 2007. 

Since at the core of the algorithm being presented there is the concept of singularities, two natural 
extensions of this work consist of applying the same ideas to compute solutions for linear ODEs of arbitrary 
order, where the equivalence can be solved exactly [14], and for second order equations under rational 
transformations, perhaps generalizing the work by M.Bronstein [8] with regards to 1F1 solutions to compute 
also 2F1 solutions. Related work is in progress. 
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